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Abstract 

We study dendritic growth numerically with a phase field model. Tip 
oscillation and regular side-branching are observed in a parameter re- 
gion where the anisotropies of the surface tension and the kinetic effect 
compete. The transition from a needle pattern to a dendritic pattern is 
conjectured to be a supercritical Hopf bifurcation. 

1 Introduction 

Crystal growth has been intensively studied as a problem of pattern formation 
far from equilibrium. |^] Dendrites are typical growth patterns found in crystal 
growth exhibited by most metal alloys and some plastic crystals. [I| The phase 
field model, which has the form of the Ginzburg-Landau equation coupled with 
a diffusion equation, is a useful model for such growth patterns.^], ^, [?], || For 
melt growth, the Ginzburg-Landau equation represents the dynamics of a phase 
transition from a liquid phase to a solid phase, and it is coupled with a heat 
conduction equation for the latent heat generated at the growing interface. 

The standard theory of crystal growth predicts that a needle crystal selected 
by the solvability condition is stable, and side branches are produced owing to 
thermal noise. || [l^] However, there are some experiments suggesting that side 
branches may appear as a result of regular tip oscillations, jy] |l3|, [TJ] Liu 
and Goldenfeld investigated the linear stability of a steady needle crystal using 
the boundary layer model and found an instability in a certain range of param- 
eter values of the surface tension and kinetic coefficients. |Q Ihle numerically 
investigated the competition between the kinetic and surface tension anisotropy 
in dendritic growth using a fully dynamical front-tracking method and found os- 
cillatory growth. ||H| However, the numerical results do not seem to exhibit clear 
oscillation, probably owing to numerical noise resulting from the discretization. 
Such numerical noise effects seem to be small in the phase field model. We have 
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numerically investigated the competition between the kinetic and surface ten- 
sion anisotropy using a phase field model and found that the growth direction 
changes from (10) to (11) as the kinetic coefficient is changed. |l7j This is a tran- 
sition from a surface tension dendrite to a kinetic dendrite. Such transitions have 
been observed in experiments of viscous fingering and NH4CI crystal, [jl^, [l9| In 
a previous work, we found that complicated patterns grow from small round 
seeds in an intermediate parameter region where the two anisotropies compete. 
However, we could not find clear limit cycle oscillation of the tip growth. In this 
work, we consider the situation in which the anisotropy parameter is changed 
slowly, in a stepwise manner, and find a transition from a steady needle pattern 
to a tip-oscillating dendrite. 

2 Tip oscillation of dendritic patterns in the phase 
field model 

The phase field model in two dimensions is written 

r(0)d t p = {p-Xu(l-p 2 )}(l-p 2 ) 

+d x {W(9) 2 d x p - W(6)W'{9)d y p} 
+d v {W(0fd yP + W{6)W'{9)d x p}, 
d t u = DV 2 u + d t p/2, (1) 

where p(x, y, t) is the order parameter, A is a coupling constant, and t(9) is the 
time constant for the order parameter, with 9 = axctaa{d y p / d x p) the angle be- 
tween the normal direction of the contour of constant p and the x axis. The solid 
and liquid phase correspond, respectively, to p = 1 and p = — 1. The normalized 
temperature is denoted by u(x,y,t) = (T — Tm)/{L/C p ) where T(x,y,t),TM,L 
and C p are, respectively, the temperature, the melting temperature, the latent 
heat and the specific heat. The diffusion constant for u is denoted by D, and 
W 2 (9) is an anisotropic diffusion constant for the order parameter. The gener- 
alized Gibbs-Thomson condition (the Wilson- Frenkel formula), 

Ui = -d {e)K-P{9)v n , (2) 

is assumed for the sharp interface model, where Ui is the normalized tempera- 
ture at the interface, do (9), k, /3 (9) and v n denote, respectively, the anisotropic 
capillary length, the interface curvature, the anisotropic kinetic coefficient and 
the normal interface velocity. Karma and Rappel derived the values of do(9) 
and j3{9) in the sharp-interface limit to the phase field model as Q 

do(9) = jj{W(e) + W"(0)}, 



I r(9) 
XJ W{9) 



W 2 (9) K + JF 



1 X 2Dt(9) 



(3) 
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where I = 2V2/3, J = 16/15, K = 0.13604 and F = \/2m2. If t{6) = W 2 {6) 
and the coupling constant A is assumed to be A = (2ID)/(K + JF), the 
anisotropic kinetic coefficient becomes zero, and the standard Gibbs-Thomson 
condition is realized. Karma and Rappel performed numerical simulations 
for this case, and obtained dendrites with tip velocities and tip shapes that 
agree well with the theory based on the solvability condition. The growth law 
vp 2 =const., where v is the steady growth velocity of the tip and p is the tip 
radius, is then satisfied. 

We performed numerical simulations of the phase field model with an anisotropic 
kinetic coefficient. We found that if the kinetic effect is included, the tip ra- 
dius does not become sharper and sharper, in contrast to the case in which 
there is no kinetic effect and the relation t>p 2 =const. is not satisfied, when the 
supercooling A is close to l.Q I n tn i s paper, we study the competition of 
the surface tension anisotropy and the kinetic coefficient anisotropy. Four-fold 
rotational symmetry is assumed for the anisotropy: W(9) = 1 + e4cos(40) 
and t(9) — W(9){1 + efccos(40)}, where e 4 and ek denote the anisotropy 
parameters of the surface tension and the kinetic coefficient. The coupling 
constant and the diffusion constant are chosen as A = l.9ID/(K + JF) and 
D = 2. This choice of parameter values corresponds to the anisotropy coeffi- 
cients d (6) oc 1 - 15e 4 cos(40) and (3(9) oc 0.05 + (e k - 0.95e 4 ) cos(46>) in the 
sharp interface limit, owing to Eq. (3). 




X 



Figure 1: Needle pattern exhibited by the phase field model with four- fold 
rotational symmetry for e 4 = 0.06, A = 0.7 and ek = 0.07. 

When ek is larger than 0.95e 4 , the (11) growth becomes favorable, owing 
to the kinetic effect. Numerical simulations were performed using the finite 
difference method with Ax = 0.4 and At = 0.015. The system size was L x x 
L y = 640 x 160. The uniform undercooling u — —A was employed as the 
initial condition. Mirror-symmetric boundary conditions were used at x — 
and y = Q,L y and the fixed boundary conditions p — —1 and u = —A were used 
at x = L x . The parameters A and e 4 were fixed as A = 0.7 and e 4 = 0.06, and 
the parameter ek for the kinetic anisotropy was changed step by step. 

The process used in our numerical simulation is as follows. If the tip of 
the dendritic pattern reaches x = 600, the growth for the parameter e k is 
stopped, and the numerical data for the order parameter and the temperature 
are saved. The order parameter and the temperature profiles in the tip region 
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Figure 2: (a) Dendritic pattern with regular side branches at = 0.06, A = 0.7 
and efc = 0.08. (b) Time evolution of the interface temperature u(y, t) for 
y = 20. (c) Amplitude |Au(y)| of the oscillation as a function of y. 



(p(x, y), u(x, y) for x > 520) are used as the initial conditions for the parameter 
efc at the next step, that is, p(x, y, 0) = p(x — 520, y), u(x, y, 0) = u(x — 520, y). 

Figure 1 displays a snapshot of a needle pattern for ejt = 0.07. This needle 
pattern is growing in the (10) direction, owing to the surface tension anisotropy. 
As the parameter is increased, the growth direction changes from the (10) 
direction to the (11) direction due to the kinetic effect. In region of intermedi- 
ate values of the parameter e^, a dendritic pattern with side branches appears. 
Figure 2(a) displays a dendritic pattern for = 0.08. Regular side branch- 
ing is seen. The interface temperature u(y, t) is evaluated for each y as the 
temperature at the position (x, y) where p(x, y) = is satisfied. For this pa- 
rameter value, the average tip velocity is approximately v — 0.265, and the 
diffusion length I = 2D/v ~ 15. The system size is sufficiently larger than 
the diffusion length. Figure 2(a) displays a snapshot pattern at the instant 
that the tip reaches the position x = 600, and the effect of the boundary at 
x = L x = 640 is not yet important. Figure 2(b) displays the interface tempera- 
ture at y = 20. The interface temperature and the curvature oscillate regularly 
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Figure 3: Amplitude of the oscillation |Au(y)| at y = 20 as a function of ek- 



at this value of the parameter. Figure 2(c) displays the amplitude of the oscil- 
lation |Au(y)| = maxw(?/,i) — min it(y,i) as a function of y, where maxit(y,i) 
and mmu(y,t) are, respectively, the maximum and minimum values in the os- 
cillatory time evolution of u(y,t). This amplitude increases monotonically, but, 
it does not seem to fit a curve such as exp(cs 1/ ' 4 ) or exp(cs), where s is the arc 
length along the interface and c is a constant, as is expected from some theories. 
It seems, rather, to fit a linear curve, cs, although we do not understand the 
reason for this. Figure 3 displays the amplitude of the oscillation |Au(y)| at 
y = 20 as a function of ek ■ As seen, there is a transition from a stationary state 
to oscillatory growth, and the amplitude increases continuously from zero; that 
is, there is no jump at the transition. The transition is therefore concluded to 
be a supercritical Hopf bifurcation. This limit cycle process can be interpreted 
as follows. The surface tension dendrite grows in the (10) direction. The side 
branches grow in the (11) direction, owing to the kinetic effect. Then, the den- 
dritic tip becomes round, the radius of curvature of the tip region increases, and 
the growth velocity in the (10) direction decreases. As the velocity is decreased, 
the kinetic effect, which is proportional to the growth velocity, becomes weak, 
and the surface tension effect become relatively large. Then, growth in the 
(10) becomes dominant, and the velocity in the (10) direction increases. As the 
growth velocity increases, the kinetic effect becomes relatively large, and the 
side branches grow in the (11) direction again. In this way, a negative feedback 
acts between the tip velocity in the (10) direction and the kinetic effect, and a 
stable limit cycle appears. 

Next, we consider the transition for kinetic dendrites. If ek < 0, e^ < 
and ek is sufficiently smaller than e±, a steady needle pattern grows in the 
(10) direction, owing to the kinetic anisotropy. As ek remaining less than 
is increased, the kinetic anisotropy becomes relatively weak, and the surface 
tension becomes dominant. Then a transition from (10) growth to (11) growth 
occurs. In the simulations, the parameters A and e± were fixed as A = 0.7 and 
e± = —0.06, and the parameter ek was changed stepwise. The system size was 



5 



y (a) 



- 

- 




( 


) 100 200 300 400 500 600 



X 

(b) 




1000 1500 2000 2500 3000 
t 



Figure 4: (a) Needle pattern for e 4 = -0.06, A = 0.7 and e k = -0.088. (b) 
Snapshot of the pattern for e 4 = —0.06, A = 0.7 and e\ — —0.085. (c) Time 
evolution of the interface temperature u(y,t) at y = 20. 



L x xL y = 640 x 120. Figure 4(a) displays a needle pattern for e k = —0.088. The 
needle pattern is marginally stable for e k ~ —0.087, and it becomes unstable for 
e/f > —0.086. Figure 4(b) displays the snapshot of a pattern for e k = —0.085. 
Side branches develop in the tip region, and the growth of the side branches 
does not saturate. Figure 4(c) displays the temporal evolution of the interface 
temperature u(y,t) at y = 20. The amplitude of the oscillation of the interface 
temperature grows but does not saturate. The side branches develop and they 
become main branches if L y is sufficiently large. The negative feedback does 
not act in this case. If the side branches grow, owing to the surface tension 
effect, the dendritic tip becomes more round and the growth velocity in the 
(10) direction decreases. Then the kinetic effect becomes further weakened, 
the surface tension effect becomes more dominant, and the growth in the (11) 
direction does not saturate. 

Elucidation of the oscillatory growth caused by the competition between 
two types of anisotropics is the main result of this paper. However, we also 
demonstrate a transition to the tip oscillation of dendritic growth, where the 
surface tension anisotropy has eight-fold rotational symmetry. fll?] , |2l|] In a pre- 
vious work, we showed the existence of tip oscillation in such a system, but 
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Figure 5: (a) Dendritic pattern for the phase field model with eight-fold sym- 
metry for e% — 0.01 and A = 0.65. (b) Amplitude |Au(y)| of the oscillation of 
the interface temperature at y — 18 as a function of es- 



have not demonstrated the transition from the needle pattern to the oscilla- 
tion. The parameter dependence of the oscillatory growth is observed when 
t{9) = W 2 (8), where W(9) = 1 + e 8 cos8#. The coupling constant A is chosen 
to be A = (21 D) / (K + JF) , so that the kinetic effect can be neglected. The 
numerical simulation was performed using the finite difference method with 
Ax = 0.3 and At — 0.01. The system size was L x x L y = 600 x 90 and the 
undercooling was fixed as A = 0.65. Figure 5(a) displays a regularly oscillating 
dendrite for es = 0.01. The interface temperature is seen to exhibit limit cycle 
oscillation. Figure 5(b) displays the amplitude of the oscillation of the interface 
temperature |Ait(y)| at y = 18 as a function of e%. A supercritical bifurcation 
occurs at e$ ~ 0.045 for this value of the undercooling. There are few real crys- 
tals that possess eight-fold rotational symmetry. However, it is not so unusual 
for the surface tension to be expressed as d(9) = l + e4COs46> + e 8 cos8#, in which 
case, the rotational symmetry is four-fold, but the crystal has a component of 
eight-fold symmetry. As the parameter e4(> 0) is increased for a fixed value of 
e 8 > 0, (10) growth becomes dominant in comparison with (11) growth. Tip 
oscillation appears for a small value of e^, but, its amplitude decreases as 
is increased. Tip oscillation is not observed for e4 > 0.016 when es is fixed to 
0.01. (The transition to non-oscillatory growth seems to be discontinuous, with 
a jump, although this is not entirely clear.) It appears that the competition 
between (10) growth and (11) growth may be important in the realization of tip 
oscillation. 
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3 Summary and discussion 



We have found a deterministic transition from a needle pattern to a regularly 
oscillating dendrite in the phase field model. This was found in surface tension 
dendrites with four-fold rotational symmetry and competitive kinetic anisotropy. 
A similar supercritical transition was found for surface tension dendrites with 
eight-fold rotational symmetry without the kinetic effect. Regularly oscillating 
surface tension dendrites with four-fold or six-fold rotational symmetry without 
the kinetic effect have not yet been found in this model. Regularly oscillat- 
ing dendrites with side branches possess mirror symmetry as a results of the 
method used in our numerical simulation. The boundary condition is not, how- 
ever, essential for the oscillation. We have confirmed that a mirror-symmetric 
dendrite grows naturally in a simulation in which a small needle pattern includ- 
ing a small antisymmetric perturbation is placed in the middle of the left side 
of the L x x 2L y box as the initial conditions. These results are not consistent 
with the linear stability analysis of Liu and Goldenfeld and the simulation of 
Ihle. They predict oscillatory growth with antisymmetric side branches for both 
surface tension dendrites and kinetic dendrites. The detailed mechanism of the 
oscillatory growth is not yet well understood, and its elucidation is left as a 
future study. 
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